Disclination-mediated thermo-optical response in nematic glass sheets 
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Nematic solids respond strongly to changes in ambient heat or light, significantly differently parallel and per- 
pendicular to the director. This phenomenon is well characterized for uniform director fields, but not for defect 
textures. We analyze the elastic ground states of a nematic glass in the membrane approximation as a function 
of temperature for some disclination defects with an eye towards reversibly inducing three-dimensional shapes 
from flat sheets of material, at the nano-scale all the way to macroscopic objects, including non-developable 
surfaces. The latter offers a new paradigm to actuation via switchable stretch in thin systems. 



PACS numbers: 46.32,+x, 46.70.De, 46.70.Hg, 61.30.Jf 

Nematic glasses (densely crosslinked networks) (l|] are 
solids with a natural state of elongation along their director, 
and contraction perpendicular, that depends on their orienta- 
tional order. Accordingly they suffer large, reversible length 
change with heating [2], illumination |3j,|4|], solvent uptake, 
pH change and any other stimulus that causes order change. 
For glasses these strains can be 2-3% and of opposite sign 
(without conserving volume) along and perpendicular to the 
director. They can then exhibit spectacular effects such as 
large bend resulting from gradients of stimuli [5] (light, sol- 
vent) or from uniform stimuli (like temperature or weakly 
absorbing light) but with a director gradient through, for in- 
stance, the section of a sheet or cantilever J2-H] (usually twist 
or splay-bend). Glasses are distinguished from elastomers 
by being so heavily cross-linked that the director field only 
changes as a result of convection due to material shape change 
from the elastic strain. Nematic glass directors do not rotate 
relative to the matrix as in nematic elastomers - they are yfl 
conventional, uniaxial elastic solids, with moduli about x 10 4 
higher than those of elastomers. We are most concerned with 
the thermally or optically inspired length changes which we 
take to be by a factor X along the director and by Ar v perpen- 
dicular to the director (where v fulfills the role of a thermal 
or optical Poisson ratio, were the deformations to be infinites- 
imal). Elastic energy cost would be associated with imposed 
changes away from these new natural shapes - we find defor- 
mation fields matching these thermo/optical changes and thus 
of zero energy in the membrane limit we take. 

Director fields can be established in the nematic liquid 
progenitor phase before crosslinking and are permanently 
recorded in the solid state achieved after linkage. In fact com- 
plex, 3-D director fields for subtle mechanical response can 
be achieved in nematic glasses ^\ via holography and sur- 
face preparation. From the point of view of device design, the 
ability to "write" an initial director field into a solid so that it 
distorts in a predictable (and reversible way) into a new shape 
with applied temperature or light would be ground-breaking. 
In particular, the ability to take flat, thin sheets of material 
and turn them into prescribed, potentially complicated, non- 
developable shapes - at nearly any length scale - is highly 
sought after. Accordingly, we consider here the elastic re- 



sponse of thin sheets of nematic glass in the membrane ap- 
proximation. We are interested in thin sheets where the direc- 
tor is uniform through the thickness. In this situation, a change 
in nematic order gives rise to in-plane stretches (or contrac- 
tions). If the initial director field is not homogeneous, then this 
change of nematic order gives rise to inhomogeneous changes 
of stretch that may or may not be compatible. A key obser- 
vation here is that thin sheets can accommodate many more 
potential inhomogeneous changes of stretch compared to bulk 
3D specimens by possibly deforming out of plane. Further, in 
thin sheets, the energy and forces associated with membrane 
mode (in-plane stretch) scales as the thickness while those 
associated with bending mode (differential in-plane stretch) 
scales as the third power of thickness. So the bending energy 
is negligible compared to the in-plane stretches. All of this 
leads to a new paradigm for actuation [Q]: If we are able to 
find a director arrangement that leads to stretches that are in- 
compatible in bulk but compatible in sheets, then the change 
of order can be used to generate out-of-plane deformations 
with a very large blocking force (one that scales with thick- 
ness). This is in contrast with small blocking forces (third 
power of thickness) associated with designs involving bend- 
ing cantilevers. In this Letter, we show that +1 disclinations 
are indeed such arrangements. 




FIG. 1: A flat nematic glass sheet with an azimuthal +1 disclination 
heats to a cone, or cools to "anticones" (n). 
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Disclination defect textures are not extensively studied, ex- 
perimentally or theoretically, in nematic solids. However it 
is known that the energy minimizing deformation associated 
locally with changes in their nematic order is not compati- 
ble in some 3-D disclinations [8]. We analyze the mechani- 
cal response and elastic ground state of the most experimen- 
tally accessible disclination defect in 2-D, that with topolog- 
ical charge + 1 . We argue that it is a powerful way to induce 
shape change, specifically through the introduction of a point 
of localised Gaussian curvature. Fig. 03 shows examples of 
initially fiat, defected sheets after heating or cooling. 

We illustrate some ways director fields with this disclina- 
tion charge could be used in concert with system geometry 
for applications, such as sharp bending, twisting or forming 
non-developable (anticlastic) surfaces and cantilevers. 

Shape response to nematic defects has been investigated 
theoretically in fluid systems where the dominant in- 

fluences include surface tension, Frank elasticity and order- 
curvature coupling. Our solids are in the opposite limit 
- elastic stresses dominate over Frank elasticity for length 
scales greater than a nematic penetration depth for solids 
Z, = ^K//j, where /j and K are the shear modulus and a Frank 
constant. ~ 10~ 10 m for a glass.) Also, in our problem 
surface energies play no role, and order resists mechanically 
induced change. 

Response closely analogous to ours has been analyzed by 
Ben Amar et al in the elasticity of botanical systems when 
anisotropic growth creates internal stresses and forces pla- 



nar systems into the third dimension 111 lL 11211 as cones or 
"e-cones" (our anti-cones). They and we deal with local- 
ized Gaussian curvature and hence localized elastic stretch. 
Such stress intensification occurs in folding and crumpling, 
see Witten's review II 1311 . However, our cones and anti-cones 
are simple, not "d-cones" fllll . and our systems naturally take 
up such shapes rather than concentrate stress in response to 
spatial crowding. Indeed, except for very weak spontaneous 
distortions, our tip extent is of the order of the thickness as 
in classical simple cones 11311 - we return elsewhere to the 
core/far-field, bend/stretch energy balance. Thus our system 
differs from the cone sources that generate crumpling. 

We consider +1 disclinations. In the membrane limit in 
which we choose to work, escape in to the third dimension 
1 14J of the director field is not possible except perhaps near 
the core, and these disclinations are true topological defects. 
For +1 defects, many different textures, though topologically 
equivalent (see Fig [2]), differ non-trivially in their mechanical 
response. We first analyze azimuthal and radial textures. 

Consider a thin sheet of nematic glass whose director field 
is azimuthal around a +1 disclination defect, as in the left of 
Figure|2j and which is flat at some reference temperature, 7b. 
As the sample is heated above 7b, the decline in nematic order 
will cause a contraction of the natural length along the ne- 
matic directors with a local elongation of the natural length 
due to Poisson effects normal to them. In a free, uniform 
glass these natural length changes would be manifested by 
actual mechanical strains so that the elastic ground state is 
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FIG. 2: +1 disclinations - azimuthal, spiral and radial textures. 



achieved. Since the chosen director field is circularly symmet- 
ric, with integral curves simply concentric circles centered on 
the defect, clearly the sheet has a problem accommodating this 
change in the natural lengths as circumferences shrink while 
the corresponding natural radii grow. 

Fortunately, in the membrane approximation where we may 
neglect bending energies, there is an obvious geometric solu- 
tion that allows the nematic glass to respond to the imposed 
thermal strain without paying the high energetic cost associ- 
ated with elastic compressions and expansions relative to this 
changed state - a cone. This may be seen intuitively by keep- 
ing track of the deformation of a circle of material centered 
on the defect. At 7b the sheet of nematic glass is flat, and 
the circle maintains the familiar perimeter-to-radius ratio of 
271. However, as the temperature rises, the perimeter wants to 
change by a factor of the thermal deformation gradient along 
the director, P — > P* = XP, where here A- < 1 . Meanwhile, 
see Fig. Q] the material (in-plane) radius is changing as well 
due to the thermal/optical Poisson effects associated with the 
perimeter's change, ro — > r 1 = X~ v tq. For glasses v is in the 
range 1/3 to 2 |3J]. Together, these transformations imply that 
upon heating, a circle of material on the sheet remains circu- 
lar, but adopts a new in-material perimeter-to-radius ratio, that 
of a circle enclosing a cone's tip (Fig.Q] 1st panel): 



P' 



2nX 1+v r' 



2k r 1 sin (j) 
sin" 1 (X 1+v ) 



(1) 



with (j) the cone opening angle. One can think of r'sincj) as 
the embedded radius. The localized Gaussian curvature asso- 
ciated with the cone tip is thus 2tt( 1 — sincj)); circles enclosing 
the tip of a non-developable cone no longer have the ratio 2% 
of the perimeter to the in-plane radius, unlike circles on the 
cone but not enclosing the tip. Away from the singularity, 
there is no Gaussian curvature and hence no shape-induced 
elastic compressions or extensions. Therefore our sheet of ne- 
matic glass responds to temperatures above 7b by deforming 
out of plane, breaking up-down symmetry in the process. The 
opening angle varies sensitively with small strains; a 2%-3% 
contraction with v <~ 2 achieved over 90°C J2] gives (j) <~ 70°- 
66°, such strains also achievable for modest illuminations in 
photo-glasses J3]. 

But what if we cool the azimuthal sample below 7b, in- 
creasing the nematic order relative to the reference state? The 
arguments relating the local changes in length along and nor- 
mal to the nematic director go through unchanged, however 
now X > 1, invalidating our previous ansatz cone solution - 
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now the perimeter is "too long" for the in-plane radius and a 
more complex deformation must result. 

We consider deformed surfaces where the height varies lin- 
early with the distance from the center or the defect. Con- 
sider a circle of radius ro centered at the defect in the fiat 
undeformed reference plane. After deformation, this cir- 
cle becomes a curve described by {r(ty),$,h(§)} in cylin- 
drical coordinates {r,§,h) in the deformed configuration and 
parametrized by the azimuthal angle (]). If the material is in the 
minimal energy state unstretched from its new natural config- 
uration, this curve has to satisfy two constraints. First, it has 
to have a constant distance R = Ar l >o from the origin (de- 
fect) and second, its length has to be equal to P = 27tA.ro = 
2tO) +v R. Thus, the curve has to lie on a sphere of radius R; 
see the trajectory in the last panel of Fig. [T] Further, since 
the length of the curve P is greater than the length 2nR of the 
great circle of the sphere, this curve has to oscillate. In other 
words, the deformed surfaces oscillate azimuthally and these 
oscillations grow linearly radially as shown in Fig.Q] We call 
them anti-cones. We also note, 
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i \c/(j) J ' \d§ 
The simplest possibilities for h are 

h(r;A,n) =Arsinn(j) 



+ r 2 



(2) 
(3) 

(4) 



for antinodal line angle a = tan -1 A (amplitude in effect; see 
Fig. [TJ and integer n (so that the curve is closed). Plugging 
this ansatz into eqn (0 gives the relationship between the r 
coordinate and (j) at constant R needed for the perimeter: 



'\ +A 2 sin 2 « 



(5) 



Returning this relation and the form of h to eqn. (O for the 
perimeter gives, grouping factors of R and simplifying, 

P(R) = 2%RI{n,A) (6) 
where 7(n,A) depends only on the scale A and the state n: 



I(n,A) 



1 



n 2 A 2 cos 2 2nnu 
1 +A 2 sin 2 27Cnt/) 2 1 +A 2 sin 2 27t«M 



Connecting the radius and perimeter as in eqn fl} gives here 
I(n,A) = X 1+V (T — 7b); as temperature and hence sponta- 
neous distortion changes, so does the character (that is, A and 
n) of the anticone. The negative Gaussian curvature localized 
at the apex of the anticone is 2jc(1 — 7). As appropriate, 7=1 
for n = - the surface is a flat plane. Otherwise 7 ranges from 
1 to \n\ for A = — > °°, see Fig. [3] 

The analogous cosine solutions simply give rise to rotated 
versions of the same surfaces for all n ^ 0, recovering the 
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FIG. 3: The behavior of l(n,A) for different values of n and as a 
function of A. When A is large, I(n,A) approaches n. 



conical solution discussed earlier for n = 0. In this case 
I(A\n — 0) ranges from 1 at A = to at A = °o, as re- 
quired. The behavior of / for n 7^ is encouraging - our trial 
solutions yield precisely the geometries that accommodate at 
zero stretch energy a cooling of our azimuthal +1 defect be- 
low 7b. However, since each individual surface is limited to a 
maximum perimeter-to-radius ratio of 2jt | n | , we would expect 
interesting transition behavior as cooling leads to strains re- 
quiring ever more crumpled geometries, with transition states 
characterized by simple Fourier combinations of the "pure" 
surfaces or, in extreme cases, exotic surfaces that are multiply 
re-entrant in (j). Figure Q] shows an n = 3 anti-cone where the 
crumples take up more perimeter than the more slowly varying 
n =2 anti-cone. See 11211 for an analysis of large amplitude 



anti-cones where the functions /i((j)),r((j)) may not be single 
valued. We conjecture that all minimal membrane energy so- 
lutions of these defects are anti-cones. 

A new paradigm for actuation now arises. Although we are 
dealing with thin sheets (and later cantilevers) in the mem- 
brane approximation where shape change has low bending 
energy, large forces can be exerted by switching on and off 
Gaussian curvature via stretch modes that arise if the natural 
and imposed geometries are in conflict. 

For radial textures, Fig. [2] the roles of the direct ther- 
mal strain and Poisson strain are swapped. Heating above 
7b now requires a shortening of the radial length due to de- 
creased nematic order while the azimuthal direction expands 
from the corresponding Poisson effect - equivalent to low- 
ering the temperature in the azimuthal texture. Hence, each 
of these textures behaves as the other under the mapping 
T - To T - T. 

Genetically a +1 disclination defect has an intermediate an- 
gle 8 of the director with respect to the radial vector from the 
defect core. As a result, the effect of the direct and Poisson 
strains are now mixed along circles and radii centered on the 
defect. Curves along which the material feels the maximal 
effect of the direct strain and none of the Poisson strain, and 
vice versa, are the integral curves of the director field and its 
normal complement, respectively. For +1 textures with < 
8 < 7t/2, the integral curves are logarithmic spirals instead of 
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FIG. 4: A +1 disclination texture with intermediate S = 45°. The 
integral curves of the director field (heavy, red) and its complement 
(light, blue) are logarithmic spirals. 




FIG. 5: Disclinations in cantilevers also cause unusual thermal or 
optical response, (a) a strip cut from a +1 texture conically deforming 
(with <]) = 45°) bends with a conical cusp. Anti-conical deformations 
are even richer (here a = 45°), with (b) anti-clastic bending (n = 
2, anti-nodes aligned with the cantilever axes) (c) pure twist (nodes 
aligned) or (d) curvature reversal (n = 3). 



symmetry breaking does not occur for the anti-cone solutions, 
as they do not break up-down symmetry. 

Consider advantageous ways to use such sheets other than 
as cone/anti-cone machines. Simplest would be to break the 
circular symmetry by cutting out a cantilever of material. 
Strips cut from these defected textures can display a wide 
range of rich behavior, Fig. [5] cusped with sharp bending, 
pure twist, curvature reversal, and combinations of these de- 
pending on the orientation of the cantilever with the defect. 
Blocking forces would depend on the dimensions of the (anti) 
conical region compared with the length of the arms. Bound- 
ary conditions are very important, and the cut out material 
must encompass the defect core. Such active response in strips 
could be used as a light activated stirrer, actuator, swimmer, or 
perhaps as thermally-sensitive, simple machines. 

"Frozen-in" disclinated director fields in responsive ne- 
matic glasses are rich and promising systems. To exploit 
them to their fullest however, will rely on understanding the 
mechanical response of all the defect charges, along with 
how interacting multiple defects influence the resultant strain- 
mediated shape change. Such understanding would allow 
blueprinting an arbitrary three-dimensional shape in a flat 
sheet and switching it on at will. 

CDM, KB and MW acknowledge support from the EPSRC. 



simple circles and radii, Fig. [4] An analysis we present else- 
where shows that the spiral angle 8 evolves with spontaneous 
deformation X. Defining b = cot 8, one finds that b — s- b/X +v . 
For instance for heating where X < 1, then the angle 8 de- 
creases - the spiral tends more to the radial direction. Also a 
mismatch between radius and perimeter as considered above 
arises: cones or anticones form depending on whether the ra- 
dius or perimeter grows relative to the other. For instance for 
an initial 8 > n/4 and hence b initially < 1, the spiral gives 
a conical response for Xe (b 2/ ^ +v \l) and anticones other- 

wise. The cone angle is q) = sin Jp^ • 

There is a new effect, however - because of the spiral an- 
gle of the director relative to the radial direction, material un- 
dergoes motion with an azimuthal component. Furthermore, 
if 8 and A, are such that an out-of-plane deformation into a 
cone is required, then this rotation, combined with the mate- 
rial's spontaneous choice to form an "upward" or "downward" 
pointing cone as the material moves from two to three dimen- 
sions, leads to a spontaneously broken chiral symmetry. This 
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